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Abstract 

We present a macroscopic model of the decay of a coherent classical scalar field into statistical fluctuations through 

-^-the process of parametric amplification. We solve the field theory (henceforth, "microscopic") model to leading 

^ order in a Large N expansion, and show that the macroscopic model gives satisfactory results for the evolution of 
(N ■ 

the field, its conjugated momentum and the energy momentum tensor of the fluctuations over many oscillations. 
Q-i 

^ The macroscopic model is substantially simpler than the microscopic one, and can be easily generalized to include 
C ■ quantum fluctuations. Although we assume here an homogeneous situation, the model is fully covariant, and can be 
^ applied in inhomogeneous cases as well. These features make this model a promising tool in exploring the physics 
of preheating. 
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Introduction 

The eras of pre and reheating after inflation (111) are generally regarded as "probably the most vio- 
lent of putative phases in cosmic history"(d). During them, the quantum (maybe effective) degree of 
freedom describing the inflaton decays into quantum and statistical fluctuations of both gravitational and 
matter fields. The main decay mechanism is parametric amplification of the fluctuations by the coherent 
oscillations of the inflaton. 

A full description of these phenomenon therefore requires an understanding of the quantum field theory 
of parametric amplification ([3]) , placed on a curved background, and includingthe backreaction on the 
geometry (lO]). Although there has been substantial progress in later years many 
; remain, such as the relevance of nonlinear effects ( fill ) and the generation of super-Hubble 



open questions 



perturbations ( 11121] ). 

Further progress is hampered by many reasons, among which we believe the most important one is 
that we do not really know the correct microscopic theory of inflation. For this reason, it is necessary to 
consider a large number of competing scenarios, giving divergent results for such a complex phenomenon 
as reheating. 

However, this variability is limited by two factors. First, for the purposes of cosmology we do not 
really need a detailed description of the process. Typically what is required is the overall evolution of the 
inflaton field on large scales, the final temperature of the radiation field - which dominates the specific heat 
of the Universe after inflation - and an understanding of the time scales involved. We may entertain the 
notion that different models may agree on these very coarse-grained observables, while diverging under 
more sophisticated probes. 

Second, while diverse, the models to be considered are not arbitrary. They must be consistent with 
such general principles as causality and the Second Law of thermodynamics. We know from relativistic 
field theory that even these very general principles put nontrivial constraints on macroscopic behavior 

(0,0]). 

These observations suggest that it may be possible to investigate the model independent, or at least 
robust, features of pre and reheating by replacing the full microscopic models by simpler, macroscopic 
models answering to the same constraints. The macroscopic model must be consistent with causality and 
the Second Law, fully covariant, respect the conservation laws in the microscopic model, and reproduce 
its equilibrium behavior. At the same time, the macroscopic model should be a consistent hydrodynamical 
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model by itself ( lllSil ). As we shall see in the following, these requirements alone essentially define the 
macroscopic model. The remaining freedom concerns the values of the parameters in the theory, which 
must be found from (numerical) experiment. 

An indirect confirmation of the feasibility of this approach is the success of a simple fluid model in 
reproducing some features of reheating (" 11161 '). The requirement of full covariance must be stressed, since 
assuming a Friedmann - Robertson - Walker background or low order perturbations thereof is inappropriate 
to analyze the evolution of super- Hubble fluctuations. 

In this paper we shall test these ideas by studying the decay of a coherent self-interacting field into 
statistical fluctuations, both from the microscopic equations of motion, and through a macroscopic model 
built according to the principles we have set up. 

In order to obtain a compelling result, we have chosen a well known problem, namely the 
decay of an 0(N) invariant scalar field, studied to leading order in a Large N approximation 
(0,lll8l,||3],Q,H],E],|23]|24]). In a certain sense, the leading order approximation is harder to 
study within our approach than higher approximations would be. The leading order system is Hamiltonian 
(see below) and therefore it is not associated to entropy production. If integrated long enough, it shows 
revivals. The macroscopic model, being a very coarse - grained version of the microscopic model, has a 
nontrivial entropy production rate, and no revivals. Therefore any agreement between them can only hold 
for times which are short with respect to the recurrence time. Higher order approximations thermalize 

n n 

(1123]) and it is possible to introduce a growing entropy already at the microscopic level (||2^), making 
them behave more closely as the macroscopic model does. 

Besides restricting ourselves to the leading order theory, we make further simplifications. We shall 
assume homogeneous initial conditions in Minkowsky space, neglecting the geometrical aspects of re- 
heating, and we shall consider only a classical field. This later approximation is sufficient for the study 
of reheating, as recently shown by more complete analysis (|27^). We assume that at the initial time the 
fluctuations are in thermal equilibrium, although not in equilibrium with the scalar field. We disregard 
initial time singularities ( jisb . We have presented elsewhere the general construction of the macroscopic 
model (1 29]). The difference between this paper and those is that here we lay the stress on showing 
that in this way it is possible to reproduce (in a much more economical way) the overall behavior of a 
prescribed microscopic model. The macroscopic model is based on Geroch's "Dissipative type theories" 
(DTT) (|Q]). 

This framework ensures covariance and causality. To determine a concrete model it is necessary to 
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specify one thermodynamic potential and the form of the entropy production. The former is found by 
requiring that the energy momentum tensor as derived from the macroscopic model reproduces the expec- 
tation value of the microscopic energy momentum tensor at the initial time. The entropy production is 
severely restricted by demanding that the macroscopic model reproduces the equation of motion for the 
mean field and its momentum at the initial time, and that a vanishing field is a stable fixed point of the 
macroscopic equations. This leaves only one decay constant undetermined, which is found by fitting to 
the numerical results. 

As we shall show, the macroscopic model succeeds in reproducing the evolution of the mean field 
and the equivalent fluctuation temperature on time scales short compared to the recurrence time, but long 
enough that initial correlations among fluctuations are washed out. This result suggests that this approach 
could be useful in exploring fully covariant, nonlinear reheating scenarios. 

The rest of the paper is organized as follows. 

Next Section is devoted to the microscopic model and its solution. We gather here some known results 
which are necessary for comparison to the macroscopic model later on. We derive the effective action 
from first principles and take the limit ^ — > to obtain the energy-momentum tensor of the field and its 
thermal fluctuations which is our quantity of prime interest. 

Section III presents the macroscopic model. After a short review of relativistic thermodynamics 
(iQl^lQ]) we present the basics of the DTT approach. We then discuss how to incorporate within these 
formalism the fluctuation field, the mean field, and their interactions. We proceed by identifying the ther- 
modynamic potential from comparison of the respective energy momentum tensors at the initial time. We 
show that the model is causal. 

In Section IV we complete the derivation of the macroscopic model by analyzing the entropy production 
rate. We present our results in the last section, where both models are compared in detail. 

A brief Appendix on DTT theory is attached. Appendix B contains some specific calculations. 

I. THE MICROSCOPIC MODEL 
A. Equations of motions 



The starting point for our microscopic model is the 0{N) invariant action (| 
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I { - 2 [2^V + } d'^x (1) 

which represent a autointeracting field under the potential 

v-(0VO = ^^W + ^a(0vO' (2) 

the /2\ term is for convenience only, producing no change in the equation of motion. It is customary 
to define a new field by adding a constraint 

[X-A((^V + ^VA)]' 



2A 

which leads to a new action 



(3) 



S2^N j + 1^ - \x^V - ^x} d'x (4) 

Representing the expectation value of the field with respect to the initial state of the theory by 

< >= 0, <x>^K (5) 
we can write the fields as the sum of mean fields and fluctuations 

if' = (l)^ + (6) 



X = K + K (7) 
Keeping only the leading term in the large N approximation we can write the effective action 

^[0^ K] = S2[<l>\ K] + t—Tr In j {d,d^ - K) 

= S^[(t>\K]+i^Trhi[-]^G-/{x-z) (8) 

Variation of the effective action gives the equation for the classical fields (jf 

- d^di^ct)' + K(f)'^0 (9) 

and for the K field, 

K^m^ + ^<l>'<f>' + ^GF{x,x) (10) 



It is helpful at this point to rotate in the internal space so that 0' = for i E [i, N — 1], and 



kN 



Albeit it is true that the longitudinal and transverse fluctuations are different, to first order in 
the Feynman propagator for the fluctuations is given by 



which is solution of the following equation: 

i-df.d" + K) Gf (x, x') = -iM (x - x') 
If we assume a homogeneous initial state, it is convenient to introduce the Fourier expansion 



(27r) 



exp {ik ■ xj 



Ukit)al + U*,it)a[ 



The normalization for the modes is 



(11) 



(12) 



(13) 



w[u;{t),Uk{t)] = -zh 

^ [/) i 5'] = fg — fg being the Wronskian. In the homogeneous case, 



K(f) = 



with the corresponding equations of motion for the mode function Uk{t) : 

d'Ukit) 



dt^ 



+ uji{t)Uk{t) = ; ujt{t) = \k\' + K{t) 



Substituting dU) en (HI]) gi 



ives 



Gf{x, x') = Gv{x, x') + Gt{x, x') 
where we recognize a temperature-independent (or vacuum) part 

d^k h 



Gvix, x') 



(14) 



(15) 



(16) 



(17) 



exp 



(27r)3 2LUfc(0) 
and a temperature -dependent part 



ik-{x-x') {Uk{t)Ul{t')Q{t~t') + Uk{t')Ul{t)Q{t' -t)} (18) 



Gt{x, x) 



(27r)3 2a;fc(0) 



{exp [^^- (f +f)] [Uk{t)Uk{t')gk + Ul{t)Ul{t')gl] 



exp 



)k ■ (f - x')] [Ukit)u;{t') + u*{t)Ukit')] 



(19) 
(20) 
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where and gk represent the initial statistical mixture: 

[ai:ai) = i2nfg,6'ik-k')6'^ (21) 



In the coincidence limit, 



(4* 4,) = {2wy'n,5''{k - k')S" (22) 



G.fe:.) = /|^^|f/.WP (23) 



and 



(27r)3 2a;fc(0) ''''^ (^'^ ' ^) Pk{t)Uk{t)gk] + 2 |t/fc(t)|'nfc} (24) 
We will assume a thermal bath of particles initially 

- ^ exp(|i)-l 

^fc = (26) 

with 

T = (27) 

the initial temperature of the bath. The lowercase-index / is explicitly used to prevent confusion later on 
and indicate the temperature of the initial bath of the fluctuations. 

1. Classical Thermal fluctuations, energy -momentum tensor and Hamiltonian formulation 

A very successful strategy to deal with reheating in the strongly nonlinear regime has been to describe 
the inflaton field as purely classical [31]; the rationale behind this approach is that there is a rapid tran- 
sition to semiclassical behavior during inflation [32] and thus, if one succeeds in finding a semiclassical 
description of fluctuations produced by the inflaton decay, the decay itself can be described via classical 
equations of motion Is^. In this paper we will focus on the particular case of classical thermal fluctua- 
tions. This is a case that is easier to implement than the quantum one while allowing us to confront the 
basic issues. We then take the limit ^ in eq.(l23t and (|24ll to obtain 

limGv{x,x) = lim / . . . |f/,,(t)|' = (28) 



?i-o ' ^ n^oj (27r)3 2cjfc(0) 
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lim Gt(x, x) 



lim 



n 



oJ (27r)3 2^fc(0)exp(^) 



-2\Uk{t)f 



(2vr)3 2ujI{0) 



(29) 
(30) 



Thus, in this classical limit, we find 



K = m^ + ^A02 + XksTf I 



(2vr)3 2c.|(0) 

The energy-momentum tensor can be computed using dH) and 



\Uk{tW 



(31) 



(32) 



specializing thereafter to Minkowski or directly by taking the expectation value of the classical energy 
momentum tensor. Either way, we find (writing only the nontrivial components) 

2^^ (2^)3 2^,(0) I 



T' 



00 



Uk{t)Ul{t) + \k?\Mt)[ 

4 

m 
2A 



(33) 
(34) 



and (no sum over i ) 

't 



1 V 



+ kBTf 



d^k 



(27r)3 2^fc(0) 



m — 



2A 



(35) 
(36) 



These expressions are the total energy density and pressure for our system. To integrate numerically 
our equations we need a finite number of variables. Since the integrands are manifestly isotropic, we 
perform first the integration over the angular variables. In problems with spherical symmetry like this 
one, this procedure leads to a better approximation of the integral than approximating the 3-D integral as 
a triple sum over Cartesian coordinates. The remaining integral is written as a finite sum: 



T 



00 



^ k=i 



m 
2A 



(37) 
(38) 
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1 > 



^ k=l 



UAk-kmik-kit) - l\kAkf \UAk.km 



m 
2A 



where we defined 



ak = 



and we set a cutoff frequency 



kmax = A — NkAk. 



(39) 
(40) 

(41) 
(42) 



Nk is the total number of modes and is the spacing. It is straightforward exercise to verify that we 
can now rewrite the whole system as a Hamiltonian system (H has units of energy density). Writing the 
complex modes Uk{t) as real and imaginary parts 



Uk{t)^Ur{t)+iUi{t) 



(43) 



and writing p = the Hamiltonian is 



1 



The equations of motions for the fluctuations are 

dH 



dUr 



2kBTfa^Ak-k)^k (l Afe • k\^ + k) Ur 



(44) 
(45) 



(46) 



and 



dH 



1 



Ur = ^ = T-^r 

oYlr 2kBTfa(Ak-k)Ak 



dH 



2kBTfa^Ak-k)^k{\^k-k\^ + K)Ui 



(47) 



(48) 



dH 



1 



Slli 2kBTfa(Ak-k)Ak 
which also give the definitions of 11^ and Ilj. The equations for the scalar field are 



(49) 



P 



(50) 
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and 

p=-K(j) (51) 

Initially, we have [/^(O) = cos(5), Ur{0) = ujk{0) sm{5) and Ui{0) = sm{5), Ui{0) = -cjfc(O) cos(5) 
where S is some random phase. K is a Lagrange multiplier and its corresponding Hamiltonian equation 
gives the (discretized) gap equation (BTT) . At t = this read 



A 

K{0) = + -</)2(0) + XknTf J2 ^k^A.-k (52) 
^ k=i 

= m' + -02(0) + ^^^^ V ^ (53) 



an equation to be solved numerically to extract K{0) as a function of the initial conditions and param- 
eters. Note that discretization and imposing a cutoff take care of both ultraviolet divergences and initial 
time singularities. No further renormalization will be needed. 

B. The energy-momentum tensor at t=0 

As we are going to see later, in order to complete the macroscopic model we need to know the exact 
form of the interaction potential. It would be very useful to have an analytical expression for the total 
energy density and pressure in the microscopic model. Our aim here is then to find an expression that is 
analytical, manageable and a good approximation of our exact microscopic model. 

Our scheme is the following. If one considers a dissipative fluid as a mixture of two simpler fluids, 
most of the information specific to each one will be lost after the mixing. We can only rely on quantities 
which can be computed using the total energy density or pressure. However, at the initial time, there is 
some specific information about each fluid that one can assumed to be known from the initial conditions. 



On one hand, we know how to describe a classical autointeracting field as a fluid (112911). On the other 
we know how to theoretically describe fluctuations (quantum or thermal) as a fluid also. We can then 
derive the contribution coming from the interaction potential, at the very least at the initial time. We will 
suppose that the functional dependence on the temperature is the same at later times. Comparison between 
the dynamics of both models will validate this hypothesis. Under these premises, the three quantities of 
interest are: 
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2A 

2 

Pit) = |- + ksTf 



dk 



m 

2A 

A;2 



,2 
k^ 



(54) 

(55) 
(56) 

(57) 
(58) 



Where H =< > represents the total energy density of the microscopic theory and P =< T^^ > is 
the total pressure. The actual temperature of the fluctuations at later times should be obtained using the gap 



equation and some suitable guess on the modes. Using the initial conditions 

^fe(O) y |C/fe(0)| = ^C/2(0) + C/2(0) = 1 we have 



2 



k^dk 



H(o)+^^} 



+^(x(0)-m^) (x(0)+m2) + ^ 



m 
2A 



fc2 



(59) 
(60) 
(61) 

(62) 
(63) 



47r2^i(0) 

where we write explicitly the functional dependence of K{t — 0) with 0o and we use the notation 
H{t — 0) = Hq and similarly with the pressure. Let's compute first K(0, 0o)- We can rewrite it as 



ik^ + K{0)) 



(F + ^(0)) 



(64) 
(65) 



That is 



dk 



(66) 



The last integral is ill-defined, the culprit being the model itself: it is the ultraviolet divergence which 
has its root in the classical equipartition theorem. In our model, the divergence is controlled by assuming 
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a finite cut-off A. Alternatively, the theory can be renormalized via renormalization of the mass ( 114 ill '). 
Integrating we found 

1 + ^^r^ / ^ arctan ^ ^ 1 K{0, 0o) = + 1a0^ + ^A (67) 

A numerical solution will yield K{0, 0o) • Note that, in the limit T ^ we have K{0, 0o) = rn^ + 1-^00- 
Let's now turn to the total initial energy density. We have 



T + ¥^^' + ¥^ -^^(0' '^o) + (^(0' '^o))''' -^^^^ ^ 
2 67r^ 47r2 \ \ . /_ft'(0, 



(68) 

indicating explicitly the dependency on the temperature. A similar calculation yields for the total initial 
pressure: 

= y + ifS^' " ifS (-^^^°' '^"^ + ^ — ^ 



K{0, 



2A 

It is possible to simplify these two expression using eq. dFTl : 



(69) 



Horn = -p + + -m 0o)0o + — X 2X- ^^^^ 



uv /; 2^ ^ 187r2 6 ^ 3A 6A ^ 

The energy density of the fluctuations Pg, is obtained by considering the above results at po = 0o = that 
is 

- "6^^ + A 2X- ' 2A ^^^^ 

- ^^^^A3-^(ir(0,0)-m2)' (73) 



67r2 2A 



and 





A^- 


2 

m 


187r2 




ksTf 


A3- 


1 


187r2 


6A 



4 

m 



(a^(0, 0) - m^) (k(0, 0) + 3m2) (75) 
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Note that pf and pj are only defined up to an additive constant. We used that small indeterminacy in 
their definitions in order to have a natural limit for pf and for T/ ^ 0. Also note that i^'(0, 0) indicates 
the solution of eq. (l67t with 0o = 0. In the case Tf 0, we find 

= f + V(M + (77) 

using © and 

Pf = = Pf {Tf = 0) (78) 



II. THE MACROSCOPIC MODEL 



A. Klein-Gordon field and fluctuations in the DTT framework 



1. Relativistic Thermodynamics and DTT 

In order to motivate our formalism, we will present a very brief summary of relativistic thermodynamics 



( OSl] '). Recall that the basic thermodynamic relations are coded in the Euler equation, which gives the 



entropy 

S=^[U + pV-fiN] (79) 

and the first law 

dU = TdS - pdV + fidN (80) 

which ensures conservation of energy. Taking differential of dTOb and subtracting dSOb we obtain the 
Gibbs-Duhem relation 

Vdp = SdT + NdiJ (81) 

from which we deduce 

S dn N dn 

(82) 



dp_ 



_ N _ dp 



We transform this into a covariant theory by adapting the following rule (|35]): 

1 . Intensive quantities (T, p, fi) are associated to scalars, which represent the value of the quantity at a 
given event, as measured by an observer at rest with respect to the fluid. 
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2. Extensive quantities (S, V, N) are associated to vector currents S"-, u"-, N'^ such that given a timelike 
hypersurface element dJ^a = riadT., then — X^c/Eq is the amount of quantity X within the volume 
dH as measured by an observer with velocity n°-. If further the quantity X is conserved, then 
X"-^ = 0. The quantity associated to volume dS^ is the fluid four-velocity, and therefore obeys 
the additional constraint = —1. 

3. Energy and momentum are combined into a single extensive quantity and associated to a tensor 
current T"^. The energy current, properly speaking, is = —T^^Ub- 

According to this rule, eq. (17^ translates into 

TS"" = -T^'^Ub + pu" - /iiV^ (83) 

which we rewrite as 

5'^ = - PhT"^ - aN" (84) 

where we introduced the affinity a = fx/T , the thermodynamic potential x° = pP"", and the inverse 
temperature vector = T^^u"". We now introduce the concept of a perfect fluid, that is a system whose 
energy-momentum tensor takes the form T"^ = g^'' + u"'u''{p + p), where p is the energy density as seen 
by an observer moving with the fluid. Usually this is not sufficient to characterize completely the fluid, 
and another equation appears in the form of a conserved current iV"^ = 0, where TV" = nu"", n being the 
corresponding density as seen by a comoving observer. 

Without friction, that is all there is to it. When dissipation enters, things become far more complicated 
and subtle. One would like to obtain an equivalent to the Navier-Stokes equations that describe adequately 
the fluid, even when the speed is near that of light. The proposals of Eckart and Landau where found to fail 
to preserve causality and have stability problems (I isl. ll^ '). To cure that defect, Geroch presented ([30]) 
a set of equations of hyperbolic character (instead of the parabolic form of Navier-Stokes) that enforces 
causality. A brief review of this formalism, is presented in appendix 1. 

In Geroch 's approach, T"^ and are assumed to be derivable from a generating function x"" '■ T"-^ = 

^ and ^'^ representing now the dynamical degrees of freedom of the theory, x"" can be further simplified 
since, as a consequence of the symmetry of T"^, we have 

X'' = ^ (85) 
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That is, all the fundamental tensors of the theory can be obtained from the generating functional x- As an 
example, a perfect fluid is obtained if x = x(C, A*), where ^ = V~C"'^a- Simple differentiation gives 



(86) 



allowing the identification 



P = Q§ (88) 

2. Thermal fluctuations as a fluid 

Thermal fluctuations could be seen as a particular case of the preceding section; the only complication 
is how to implement the cutoff. In the case of the fluctuations, there is no chemical potential, and we are left 
with /3" as the thermodinamical degree of freedom. However, the cut-off condition selects a special frame. 

Let T" be the 4-vector associated to this frame. We may think of T° as the 4-velocity of the observers in 
whose rest-frame the cut-off is imposed. Our generating potential will take the form Xf = Xf (C) ^) where 



^ = V^^^ and = -pa^"- Then 



X" = 41 - T"^ (89) 



and 



( -7^ 1 + ( 7^ 1 (90) 

(91) 



^ ^ ^ d^dv dv 

For our goal, it is sufficient to choose the following ansatz for x 



Then 



Therefore 



xi^,v) = XiiO + X2{v) (92) 
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but 

pf^W^^^^ ^^^^ 

To solve explicitly, we identify { with j3. Then 



Xi = -j PPfdP (96) 

and 

+ J ^^^^ 

As an example, if we used only the leading term we have 

A3 A3 
= 18^ ^ = -18^^ ^^^^ 



and, using v — —PaT"' — I3T^, 

A^ V 9^X2 



(99) 



67r2/3 (5 dv^ 
we obtain 

A3 

X2^-^{vlnv-v) (100) 

Equivalently, one can write v — P to perform the computation, replacing /3 — > v in the expression for 
X2 afterward. 



3. DTT formulation of the the Klein-Gordon theory 

Casting an autointeracting scalar field in the framework of DTT is akin to the description of Landau of 
a superfluid as a mixture of a normal fluid and a superfluid. The super fluid has no entropy by itself and 
therefore cannot have a temperature defined for it. The scalar field, as a coherent solution of its classical 
equation of motion, will have zero entropy also and one thus expect difficulties in trying to describe it as 
a classical fluid. However, one can define a larger theory via the following equations 

= (101) 

and 

adopting the following constitutive relation for the energy-momentum tensor 




(103) 



The scalar field is introduced by writing the functional relationship between R and T parametrically as 
R = V (0) and T = V{(f)). If one forces the constraint 

jb,c-jc,b = (104) 

then = (f)'°- and we fall back on the usual Klein-Gordon theory. One can introduce new variables ^ 
and P^^s^ and a generating functional = P{c)P- (L29]) we show that, the Klein-Gordon theory can be 
reobtained if 

Xe = -^e'lnAe) + ^V(0)/3j) (105) 
with the current given by = -^P(^c) > ^^'^ 

e = -Pic)^,a. (106) 

The Klein-Gordon equation is given by = V'{(f)), V{(f)) being the potential. Direct computation yields 
S"- = and S'^^ = as it should be. Note that in the 'rest' frame where (3'^^-^ = (/5(c), o) we see that 
the canonical momentum p = </> is given by j9 = — ^/ /3(c) . In Klein-Gordon theory, only this ratio has 
meaning. In order to break this indeterminacy, we must look at a larger framework where the Klein- 
Gordon field interacts with another fluid. In this larger context, eg. (11061) will be taken as the definition 
of the scalar field ^. The theory describing the scalar field and a perfect fluid (f-fluid) together can be 
obtained enlarging the generating potential by the addition of an interaction functional S, as we explain in 
the following section. 



B. The mixture of the two perfect fluids 

After coupling the scalar field and fluctuations, the energy-momentum tensor of field and fluid will not 
be individually conserved and furthermore the Klein-Gordon equation will also deviate from its original 
form. The total set of equations governing our theory can be written 

j^^a = V'{<P) + A (107) 
T^^;b = I" (108) 
T^fy,b = -I"" (109) 

(il)4>,a = -e (110) 

where a subindex c (/) indicates a quantity belonging to the scalar field (thermal fluctuations). j° and 
T("^ are the current and energy-momentum tensor for the scalar field. In this paper, we will work under 
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the assumption of homogeneity, which will reduce the number of equations for the fluid to four. The 
total system, consisting in the union of the two (perfect) fluids, differs in two ways from the initial q and 
c-fluids. First, the generating functional of the theory is generated not only by the sum of the generating 
functional of each fluid. There is also a third potential to include the interaction between the two fluids 

= + + (111) 

Thus 



dP{c)h 



T^:) = -^n — - (112) 



9l3if)b 

Moreover, the mixing will also generate entropy. The entropy production can be computed using the 



T^^^ = J-A11. i (113) 



equations of motion: 

VaS" = -^^-3,1" (114) 

The pressure is by definition given by 

I{ = T' (115) 

In the homogeneous case, the only nontrivial independent equation for the energy-momentum tensor is 
the a = 6 = one. It is convenient to work with 



= Tg+T(f) = ^ (116) 



-lafe rpab I rpab 
"^ab n^ab rjiab 



T^" = T(«^ - Tf/) = 2^ (117) 

where the new variables 5" = — P^j^^ and 2/? = P^^^ + To compute (II 161) and dl 17l) we need to 
propose some specific form for the interaction functional S. We use as model the following interaction 
potential 

E = F{u)v + G{u)w^ (118) 

which can be understood as a lowest order Taylor development in B'^ assuming moreover that it is depend- 
ing only on the temperatures. We have introduced the following variables: 

u = -I3al3^ = l3^ (119) 
V = -BaB"" = (120) 
w = -BaP'' = (3B (121) 
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which represent the scalars that can be made from the mean inverse temperature 13 and 
temperature difference B. Starting from (11181) we computed elsewhere ( ll29ll ) 



\ du du 



\ du^ ctM^ J \du du ^ 

dG 

+4 (p^B^ + B"!]^) —w + 2B''B^G 



It is straightforward from there to calculate 

92° 



dl3o 



'4/3^^ + 2^ + 4/3^^ + 10/3^^ + 2g] B^ 
du^ du dw^ du 



dl3i \du du j 

It is easier to work directly with the /3— derivative. We have 

dF dy/udF 1 dF 



du 



du d(3 2(3 d(3 



Therefore 



and 



d^F _ _ 1 dF 1 d'^F 
~d^ ~ ~4fP'df3^4f3^ d(3^ 



d-EP (d^F ^^d^G ^AG ^ „2 



dE^^_(ldF dG\ 



Instead of working with F and G, we will define the dimensionless / and g as follows 



F = p-'f 



That is 



G = p-'g 



d'' ff + g' 



dPo dp^ \ p^ 



B' 
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Also 

We have to find the functional form of / and g. This will be achieved by comparing with the micro- 
scopic model to which we will now turn. 

C. Computing the interacting potential from the numerical data 

After this rather long interlude, we are now ready to return to the initial quantities that describe our 
field at t = 0. Working from the approximate form of the gap equation, the total energy and the pressure, 
we will obtain an appropriate form for the interacting potential, which means that we will work out the 
specific form of the functions F and G in eq. dl 18I ). Initially, the total energy density and pressure for the 
fluid,in the microscopic model, are given by: 

rr = f + l'('#'o) + W + ^ (135) 

T: = f-V^(^„)+p, + g (136) 
To connect with the microscopic model, we will demand, that 

= Ho (137) 
V = Po (138) 

with Hq and Pq define in dTOb and (ItTT) respectively. This is a two fluid model. Each fluid possesses 
a temperature, or inverse-temperature P = [ksT)^^. We thus have here two temperatures and Tf 
representing the temperature of each fluid. It is convenient to work with 13 = (/3c + /5/)/2 and B = Pc — Pf 
. Defining 

2 

Dh = Ho-^~V{4>o)-pf{0) (139) 

2 

Dp = P,-^ + V{<Po)-Pf{0) (140) 

where p/(0) and Pf{0) are defined in eq. (17^1 and (1751 . As (po 0, numerical simulations indicates 
that the D's goes to zero as 0o- The functions Dh and Dp can be approximated with great accuracy as a 
low order polynomial (typically of order two or three). 
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imax 



Dh 



(141) 



i=0 



and 



imax 



Dp 



(142) 



Both fit were computed to the same order. The coefficients q and di can be obtained numerically using 
a least square fit. Note, for fixed (/)o, we have as T ^ 0. As an example, we show -Dj/ at fixed 

00 as a function of T (see figure ([T])). The fit is shown superposed to and is practically indistinguishable 
from the exact expression. We also show the exact Dh and Dp at fixed T as a function of 00 compared 
with out quadratic ansatz in figures and Q. The full line shows Dh and Dp defined in Eqs. (11391) and 
(11401 ). The dotted lines are the polynomials approximations (11411) and (11421) . respectively. 

1. Relating Bq to 0o 

In the microscopic model, there is only one temperature Tf. However, the description of the 2-fluids 
needs two. At equilibrium, we know that there is only one (uniform) temperature present. That is, at 
equilibrium, the difference between the temperature of the fluctuations and of the Klein-Gordon fluid 
should be zero. In this limit, the interaction between the two should go (quadratically) to zero. Turning 
now to the microscopic formulation, we saw that, as the initial values of the scalar field goes to zero, the 
interaction term goes to zero quadratically. The initial value of the temperature difference is fixed by the 
condition that the total energy-momentum tensor of the macroscopic model matches the same quantity 
of the microscopic model. As the latter depends on (pQ, this suggests that 0o and the initial temperature 
difference are linearly related 



This equations defines the 'field' temperature Tco at t = 0. Following our philosophy of simplicity, we 
choose to model the proportionality term via a power law with a as parameter, a will be chosen in 
order to ensure (local) causality. The cutoff, in inverse wavelength space, has the units of temperature in 
natural units. Moreover, it makes sense to have the interaction potential proportional to the cutoff since 
the fluctuations are proportional to the cutoff in the microscopic model (recall Eq. (l67b and the paragraph 
before it). In practice, this means that we will consider the following relation: 




(143) 
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{kBT,-kBTff (144) 

f) IF 

In the last line we evaluated Pc, Pf zero order, that is, we take Pc, Pf — P since we can only seek a 
second order term. We thus obtain the following equations: 

iTThClX 



i=0 
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250 



exact 

H 




*0 

Figure 2: Exact Dh ( eq. ((HD) as a function oi 4> utT = Tf = 14. The full line is Dh in Eq. (fT39l . The dotted 
line represent the quadratic ansatz given by Eq. (I141t . 

1=0 

with the Q and di being constant (but still dependent of the parameters of the theory, namely m and 
A). Observe that by adopting this definition we model well our interacting term in a certain range of 
temperature around Ty but we cannot now take the limit Tf 0. It is convenient to introduce the 
following definitions 

imax 

f^{T) =A''J2 CiT^'^^-"^ (148) 

i=0 
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imax 

f\T) = A" 5] diT^'^^-"^ (149) 

1=0 

We are now able to compute the interaction potential that we are going to use in the macroscopic model. 
However we must pay a price for the crossover from micro to macro. Indeed, depending on the value of a, 
the limit T — > could be ill defined. If we return to the microscopic model, our limit — implies that 
we do not expect our model to be valid at low temperature. Indeed, we do not expect to obtain valid result 
when ksTf <^ ^^^(0) = ^k'^ + K (0, 0o) (see Eq. (l25t ). This indicates that we could obtain a minimum 
temperature comparing the thermal energy ksTj to the effective mass at 0o = 0. To find T^m we consider 
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2 4 6 8 10 12 14 16 18 20 

T 

Figure 4: ^yK{0, 0) as a function of T . The intersection with the diagonal dotted line determine Tmin ■ The value 
of Tmin depends uniquely on the parameters m and A. With m = 1.1 and A = 15 we obtain Tmin — 10.198 

\J K (0, 0) versus T (see figure). For value of T less then \Jk{0, 0) then we expect that our ansatz could 
go awry. We will then define Tmin as the (unique) solution to the equation 

{0,0, Tmin) - Tmin = (150) 

with K{0, 0) the solution of eq.dSTt with (f)o = (see fig.©). With m = 1.1 and A = 15 we obtained 
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D. Computation of F and G 



The interaction potential (II 181) is defined by two functions of the mean temperature F and G, or equiv- 
alently the dimensionless / and g defined in (11311) and (I132I) . Now that we have extracted what should 
be the interacting part from the microscopic model at time t = 0, and made the correspondence with the 
variables (that is T and B) of the macroscopic model, we are ready to compute explicitly the interaction 
potential. Using (11331) and (11341) together with (11481) and (11491) we obtain: 

(f + g\ 



P (151) 



d(3^ V 



with = and = f\j3-^) as given by (fT48l) and (fT49l) . To ease the notation we will drop 

the tilde from now on. It is not complicated to turn around these equations to obtain expressions for / and 
(yf . As we explained earlier, we have to set our limit of integration in the range where we believe our ansatz 
for and t** to be acceptable. Defining (S^ax = T^j^ we obtain: 

where a and h are integration constants, h in fact will turn to be irrelevant. We have 



f + g = j3^ dx dy + ap^ + h(3^ (153) 

J a Jx 



= ^ (2/ + 3^?) - \P' dx dy t'\y) + + P% (154) 



leading to 

1 „ .. 1 „ P0max 1 r 

9 = ^P'f^ - dy f'{y) + -P'a (157) 

Replacing in the expression for / + we immediately find 

f = P'j^ dx dy t^^^y) - -P'r + -P' dy t'\y) + -ap' + bp^ (158) 



That is, using (fmi) and Sim 

pPmax fPmax In- 1 f Pmax \ 

F = J^ dx dy t'\y) - -P'r + -P dy t'\y) + -aP + h (159) 
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and 

The causality conditions involve derivatives of F and G with respect of their argument u = P"^ = T~^. 
If we denote by a prime a derivative with respect to m , we have the following relation: 

^ = (161) 

du 2 dT ^ ' 

Reverting to our old notation (that is = t^^{T) and similarly with t") we thus obtain 

F' = -\t dz^^t'\z) - \nT) - \e\T) + \§^t%T) + \Ta (162) 
G = - f\T) - - r dz\t^\z) + -Ta (163) 

2 2 JTmin Z 2 

G" = -^r='^r(T) + r \f\z) dy + -TH'^T) - -aT' (164) 

with the prime denoting differentiation with respect to u = T^^ and a lower integration limit T^m > 
was put in to ensure convergence (this problem is a consequence of the ansatz Eq. (11451) as we anticipated 
in the previous section). 



E. Restrictions imposed by causality and tiie second law of thermodynamics 

Now we can analyze the restrictions placed on the model by causality considerations. Let us recall the 
(local) causality conditions 

dF 

— > (165) 

du 

rfG /rfF^ \ <0 (166) 
du (3'^ \du I 



'dG 1 fdF \ \ dF 



^ was written before (11621) . Note that the last two conditions imply the first. Now 

du ~^ P"^ \du ~^ ) 2 Irmin z"^ ^ ^ ~^ 2 ^ ^ 
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It is simple algebra to write the last condition but the expression is very messy and not very illuminating. 
Note that the second condition is verified whenever t^^ > which is the case in our microscopic system. 
Note that b does not appear and is therefore irrelevant. The conditions are fulfilled for a wide range of 
value of a. However, we will see later that the best fit for the model happens with the value of a that 
almost saturates the conditions, that is that make the conditions almost zero for T = Tj . The value of a 
that allows the best fit with the microscopic theory was seen to correspond to a = 0. 

The first conditions (11651) is shown in figure ©, the second condition (11661) in figure (0) and the third 
(11671) in figure d?]). All are plotted as a function of T, with the choice a = 25.9. (with m = 1.1, 0o = 2*m, 
Po = 0, A = 40, A = 15 and T„i„ =~ 10.198. as given by figure dH)). 

Usually the conditions are fulfilled within a range of value for a that depends on the value of T^m.- 
ais greater then ~ 6 with Tmm — but the minimum value for a grows with increasing Tmm. Therefore, 
within a finite range of temperature ( from approximately T = 10.6 to T = T/ = 14 in the case of figures 
(13) to (171)), the three conditions are met thus ensuring causality and stability. 

ni. COMPARISON BETWEEN THE MICROSCOPIC AND THE MACROSCOPIC MODEL 
A. Matching the dynamics 

To compare the two models, we have to choose parameters in a consistent way to ensure that we are 
comparing effectively the same physical system. Some parameters are trivially chosen since we assume 
that the mass of the particle m and the coupling constant A are exactly the same in both models. We also 
have already taken for granted that the temperature fc^T/ in the microscopic model is equal to PJq = 
{Po — . The quantities that should be equal, at least initially are the total energy density, the total 

pressure and K{0, 0o)- Note that the initial condition for the field and its conjugate momentum should 
also be the same in both theories. Our microscopic system conserves both the total energy density and the 
entropy. However, since the macroscopic system is a coarse grained version of the full system, one could 
expect some generation of entropy in the macroscopic model. To describe this, we will allow a nonzero 
entropy production. The entropy current is defined as follows 

5" = - - - (169) 
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Figure 5: The first causality's condition, eq. (11651) . as a function of temperature. The allowed range of temperature 
is given when the function is positive. The value of a = 25.9, with Tmin = 10.198 and the integration constant a 
(see eq. (I153t )was set to zero. Other parameters were m = 1.1, A = 15, the cutoff A = 40 and the initial conditions 
(1)0 = and po = (see eq. (fHSl and (fl?^ ). 



As we said earlier, the entropy production is computed to be 

VaS" = -Bar - 

where the equations of motion (I107I) . (11081) and (11091 ) have been used. Defining 

A = A^ + C/3aB'' 
29 



(170) 

(171) 
(172) 







condition-&-5 



-1 




10 11 12 13 14 15 16 17 18 19 20 



T 

Figure 6: The second causality's condition, eq. (11661) . as a function of temperature. See the caption of fig. ^ for 
details. 

We can rewrite (11701) in the following manner 

VaS" = -iA^ + CI3aB'')i - [d^/S" + EB'' + FptB'/S") Ba (173) 

If one demands Va^"" = then A = E = F = (C + D) = 0. Relaxing this restrictions to allow some 
entropy production will introduce more terms in the equations of motion. More specifically, we use the 
following ansatz for the A and terms: 

A = -A^-C(3B (174) 
/o = -C^f3 (175) 
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Figure 7: The third causahty's condition, eq. (I167t . See the caption of fig. ^ for details. 

That is, 

A = -C8B - 4p + ifi) (176) 

where p is the momentum conjugate to and C has units of and A of T^. We choose to model this 
factor as follows 



c = CqTV 

A = Ao0' 
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(178) 
(179) 



Recall that (j) has units of temperature in natural units, that we adopt in this work. The (j) dependence 
is to enforce the fixed point 0c = 0. Aq is the entropy production term. Entropy production will be 
linked with the decay of the scalar field since, as a consequence of the conservation of energy, this decay 
will imply a growth of the energy density of the fluctuations. Without it, the amplitude of the field stays 
constant in time. The value of Co is chosen by demanding that the initial value of the left-hand side of the 
equation for p is the same for both theories. We thus obtain: 

We can now return to the equations (11071) to (II 101) and write them as a function of our variables. In order 



Co 



K{0, 0o) - - ^ ^^^^ Q_ ^ 1 



to integrate them numerically, it is convenient to rewrite them in the following form: 

= P (181) 

p = -V'{(f))-A (182) 

T = -Aph{T,B)-2I^h{T,B) (183) 

B = Apfs{T,B) + 2l'UT) (184) 

where A and are defined en eq. (11761) and (11771) . The functions fi{P) are quite complicated and are 
obtained upon substitution of the interacting potential in the equations of motions. Reordering to display 
the equations in a form suitable for numerical integrations, 

f, = (185) 

gigi - 9293 

f2 = (186) 



(187) 



/4 = — (188) 
where 

g, = ^Tfp'^ + T[{T)B' (189) 
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g2 = lTfp'^ + 2n{T)B (190) 



93 = -^T]p'^ + a[{T)B (191) 

g, = -^T}p'^ + a^iT) (192) 
with means that the derivative is taken with respect to Tf 

imax 

ri(T) = A" ^ Qr+^-° (193) 

i=0 

0-1 (T) =4 r \T,{z)dz (194) 

with the Q obtained by least square-fitting as explained in section 2.5.1. The integration constant a that 
appears in eg. (11531 ) has been set to zero, as discussed above. 



B. Results 



We are interested in describing the scalar field and also to be able to predict the temperature of the 
fluctuations. We do not expect our coarse-grained model to match exactly the microscopic one but at least 
we should be able to follow the decay of the inflaton and the temperature (which means that the density 
of the fluctuations should be accurately predicted). Of the various parameters that appears in our theory, 
almost all are fixed by the initials conditions or the parameters of the theory. T^m is fixed by Eq. ( 11501 ) 
and depends only on m and A. Cq by Eq. dlSOl ). Finally the constant of integration a was set to zero. 
There are two parameters that can be adjusted: Aq (see eq. ( 11791 )) and a. In both cases, these parameters 
are not fixed by the initial conditions but are chosen to match some behavior of the scalar field, namely, to 
emulate the decay of the inflaton field in the case of Aq and to match the phase in the case of a. However, 
once fixed, they cannot be changed when another set of initial values is chosen. In the examples shown 
below, Aq = 0.9. and a = 25.9. Moreover, as we have seen before, the value of a is the one that almost 
saturates the causality conditions. 

We choose the following set of parameters for the numerical simulations that we present in the article: 
m = 1.1, A = 15 (see eq. Tf = 14 (see eqdZTt). The high value of A is to ensure that the nonlinear 
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effects are important. The cutoff A = 40 (see (H^ ) with a spacing A^. = .05 (see (IISI ) which gives a very 
good approximation of the integrals by Riemannian sums. This implies that we work with 400 complex 
modes in the microscopical model. This was integrated in Fortran using the Bulirsch-Stoer method as 
implemented in the Numerical Recipes A simple Runge-Kutta adaptive step algorithm was used to 
compute the equations system for the macroscopic model. To compare we show the homogeneous mode 
of the scalar field of the micro version against the scalar field in the fluid model. 

Our objective is to demonstrate how the macroscopic model captures the essential features of the mi- 
croscopic one. In order to do so we display a number of graphics that compare the scalar field and its time 
derivative in the two models. The initial values are 0o = / * m and po = 0. For the first run, we set 
/ = 1 and the scalar field is given as a function of time. In this case, the scalar field by itself is not very 
different from the one given by the Klein-Gordon equation with rnljrj = K{0, 0). The parameter a = 25.9 
( see (I143I) '). the value that almost saturated the the causality conditions (see (I165I) . (I166I) and (11671) ') , with 
Tmin = 10.198, eq (11501 ) . In figure we show the energy density. Note that in both cases, this is 
conserved up to numerical accuracy, as it should be. The slight difference between the total energy density 
in the two models is generated by the error inherent to the least square fit approximation for pf (almost all 
the energy density came from the fluctuations). 

In figure the scalar field is depicted as a function of time for both models. For these values, they are 
almost equal up to numerical accuracy. To ease the comparison, the two graphs are shown superimposed 
in figure (flOb . In figure (fTTT) . the energy density of the fluctuations for both models is also shown. We 
recall that 

p---(t) = 'g%,T;(t) (195) 
where the are obtained by least square fit of eq. (1731 and 

r/ = (i-iB)" (1%) 

and 

pr™(^) = - Ip'it) - Vi<l>ii)) (197) 

At t = 0, the two should be equal since the interaction energy is null initially. Note that p/ is a function of 
Tf. To make the reading easier, the output in both cases was processed with a high frequency filter to get 
rid of the high frequency and to show the general trend. In the microscopic model, we define the energy 
density subtracting from the total energy density the energy density of the scalar field. 
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In the next set of figures (fig. lfT^ to (fT5l). we set / = 2 in order to obtain a more pronounced damping 
of the amplitude of the homogeneous mode of the field. Note in figure 7 that the macroscopic model 
captures very effectively the damping of the amplitude for large time. None of the parameters of the fluid 
model was changed, that is the same value of a that was used in the case / = 1 was retained, namely 
OL = 25.9, a = 0, Train = 10.198 and Aq = 0.9 in order to show the predictive power of the model.. Under 
the rather crude assumptions that lead to our model, the result is quite satisfactory. 

In figure (fT?b . we compare the time derivative of the scalar and the macroscopic model. Finally, we 
depicted in figure (fTSt . the energy density for the fluctuations in both models in the case / = 2. 

We conclude that the macroscopic model is able to describe the evolution of the field and the energy 
density of the fluctuations with the accuracy required for cosmological applications. 



Appendix A: Divergence Type Theories 

Following Geroch divergence type theories are usually described in terms of some tensorial quan- 
tities that obey conservation equations 

T"^^ = (198) 
N"^^ = (199) 

^abc^^ = Jbc (200) 

This is a simple and slight generalization of relativistic fluid theories proposed initially by Liu, MuUer 
and Ruggeri I^Sj]. In this setting, T"^ is the energy-momentum tensor and A^" is the particle current. 
Their corresponding equation simply expresses conservation of energy, momentum and mass. The third 
equation will describe the dissipative part. The energy -momentum tensor is symmetric and A"-'"^ = A"'^''; 
v4"^ 6 = and 1°:^ = 0. The entropy current is enlarged to read 

S'' = x''- ibT"^ - iN" - ibcA"^^ (201) 

I — I 

The the dynamical degrees of freedom. The following relations hold L37.1 

iV^ = ^ (202) 

j.ab ^ ^ (203) 
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Figure 8: Energy density as a function of time for the case (l)o = m. The soUd Une is the macroscopic model and the 
dash line correspond to the microscopic model. The slight difference came from the error in the least-square when 
we compute the energy density of the fluctuations,/?/, and the i nteraction energy, t^^B^ 



j^ahc ^ ^ ^204) 
Symmetry of the energy-momentum tensor implies that 

= ^ (205) 
That is all the fundamental tensors of the theory can be obtained from the generating functional x- The 
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Figure 9: The scalar field shown as a function of time. The upper plot correspond to the microscopic theory and the 
lower plot to the macroscopic theory. The initial value for (pQ = m, the parameter Aq = 0.9, the parameter a = 
and a = 25.9. 

entropy production is given by 

•5 ;„ = -I'Xbc (206) 

Positive entropy production is ensured by demanding that I^^ = M^^'^'>^'^^'> ^de, where M is negative 
definite. 

Ideal fluids are an important if somewhat trivial example. To obtain ideal hydrodynamics within the 
DTT framework, consider a generating functional Xv = Xv (^^ lA where /i = \/—^^. It is a simple matter 
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Figure 10: The scalar field of fig.® superimposed. The solid line represents the macroscopic case and the dotted 
Une the microscopic theory. 



to obtain 



edXp 



T. 



ab 



+ 



+ 



ji djj. ' jj? 

A simple comparison with the perfect fluid form of the energy-momentum tensor T"-^ 
u°'v!' [p + p] implies the following identification 

1 dXp 



P 



(207) 



(208) 



gabp _|_ 



(209) 
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Figure 1 1 : Mean value of the energy density of the fluctuations. The solid line represent the macroscopic theory 
given by eq. ( 11951 ) and the dotted line the microscopic theory given by eq. ( I197t . 



Note that the conserved current can be quite generally written as 



P = ^ (210) 



A less trivial but important example both historically and conceptually is the Eckart theory which can 
be obtained from 

XE = Xp + ICabu'^u' (212) 
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Figure 12: The scalar field depicted as a function of time in the case cpQ = 2m. The upper graphics represents the 
macroscopic case and the lower one the microscopic case. All parameters where kept to the same value as in fig. ^ 
namely Aq = 0.9, a = and a = 25.9. 

Performing a Legendre transform to the new variables ^ one obtains a system of first order 

differential equations of the form 



<B;a = I' 



(213) 

where stand for the entire collection of variables , ^ab) and similarly = (o, 0, /"^) represent 
the dissipative source; the index A thus covers 14 dimensions in our example. This first order system of 
differential equations is symmetric since 
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Figure 13: Same as fig. (fT2l but shown together. The solid line is the macroscopic theory and the dotted line is the 
microscopic theory. 



Note that we have a system of the form 



A'vi + Bv = 



(215) 



where i is a space-time index, the A* and B are k x k matrices and v is a. k -vector. Now this (first 
order) system is hyperbolic if all its eigenvalues are real; each of these eigenvalues represent the velocity 
of propagation of some small disturbance in space. These in turn propagate along^ hypersurfaces called 
characteristics whose existence is insured by the existence of k real eigenvalues [39] [40]. If the matrices 

and B are symmetric then it suffices that some combination A^Vi be definite (negative-definite given 
our choice of the signature for the metric) to insure that all the eigenvalues are real (but some could be 
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Figure 14: The same initials conditions as in fig. (fT2l but the time derivative of the scalar field is shown as a function 
of time. Solid line is the macroscopic theory and the dotted line is the microscopic theory. 

degenerate). An usual case happens when this combination reduces to A° , the vector v being the time- 
like vector (l, if ^ . In a relativistic theory one would expect hyperbolicity to be invariant under (proper) 
Lorentz transformations; in this case we say the system is causal. In our context, one would thus say that 
the system is hyperbolic if 



'^^ (216) 

is negative-definite for some temporal vector Va and the theory will be causal if this stay true for any 
temporal vector Va- 
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Figure 15: The mean value of the energy density of the fluctuations in the case (/)o = 2m. and parameters as in fig. 
(fT2t . Solid line represent the macroscopic case and the dotted line represent the microscopic theory. 

Appendix B: Computation of the expectation value of the energy-momentum tensor: 



We begin with 



2 bV 
2 bS2 



+ 



2 bV 



f 



9 ^9au 5fl 



(217) 
(218) 



using the well known identities 
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This leads to 



(220) 



(221) 



with the prime denoting differentiation with respect to x' and where 



GT(x,x')^kBTf 



d^k 1 



exp 



ik ■ (f - f )1 [Uk{t)Ul{t') + Ul{t)Uk{t')\ (222) 



n (27r)3 2a;K0) 

with the limit ^ ^ already taken care of. Let us consider first the case ii — u — O. We have 

m - Iv'' i-m + m gt(x,x')i_, = {^-&a + \d[d^ gt{x,xx=.' (223) 



(27r)3 2u;^(0) 



Therefore 



(fk 



In the case n — u — iwe have 
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(224) 
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/r/3i'1(' 2 1 -> ^ 



Thus 



/• (i3fc 1 r 
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U,{t)\"-[^K+^-\k\']\U,{t)\' 



The classical energy-momentum tensor of a scalar field is given by 

2 (5^2 



rpfXV 
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with 
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where we factored out a factor of A^. We obtained 
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2A 
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that is 
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which leads to 
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Now 
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which, using (l3Tt reduces to: 
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Therefore 
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leading immediately to the final answer 
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